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1. Introduction 

In this paper, we consider the bonndary valne problem L = L{i, U, V) for the eqnation 

iY :=¥"+ {p^I+ 2ipQi{x) + Qo{x))Y = 0, a;e(0,7r), (1) 

with the bonndary conditions 

U{Y) := W(0) + (tph + ho)Y{0) = 0, 

V{Y) := Y'{7t) + {tpH, + Ho)Y{7r) = 0. ^ ^ 

Here Y{x) = [yk{x)]i.^Ym ^ colnmn vector, p is the spectral parameter, / is the m x m nnit 
matrix, Qs(x) = [Qs,jk(x)]j^k=pm m x m matrices with entries Qs,jk{x) G iy^^[0,7r], s = 0,1 
hs = [hs,jk]j,k=T:^, = [Hs,jk]j,k=T:^, where h^jk, H^jk are complex nnmbers. 

Differential eqnations with nonlinear dependence on the spectral parameter, or with so-called 
“energy-dependent” coefficients, freqnently appear in mathematics and applications (see m 
and references therein). In this work, we stndy the inverse problem of spectral analysis for the 
matrix pencil L, which consists in recovering the coefficients in ([T]) and ([2]) from the spectral 
characteristics of L. In the scalar case (m = 1), snch problems were stndied in papers |BHS]. 

If Qi{x) = 0, then L becomes a matrix Stnrm-Lionville operator. The inverse problems the¬ 
ory for these operators was developed in works [IDHI5]. We have also stndied inverse problems 
for qnadratic pencils of matrix Stnrm-Lionville operators in [ISIE]. The case of the pencil on 
a half-line [16] appeared to be easier for investigation since this pencil has only a finite nnmber 
of eigenvalnes. In [HI, the case of a finite interval was considered and the nniqneness theorem 
for the inverse problem was proved. In the present paper, we continne this research. Onr goal 
is to provide a constrnctive procednre for recovering of the pencil L from its spectral data. 

We assnme that the coefficients of the pencil satisfy the conditions: 

(I) det(/ ± hi) ^ 0, det(/ ± H^) ^ 0; 

(II) Q|(x) = -Qi{x), QI{x) = Qo(x) for a.e. x e [0,7r], hi = -hi, ihj = -ihi, hj, = ho, 
Hq = Hq (the symbol f denotes the conjngate transpose). 

Condition (I) exclndes problems of Regge type (see [IB]) from the consideration, as they 
reqnire a separate investigation. 
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Let be the matrix solution of equation ([T]) under the conditions 

[/($) = /, ld($) = 0. We call <h(x, p) the Weyl solution for L. Put M(p) := <h(0,p). The 
matrix M(p) = [Mjk{p)]j k=pm called the Weyl matrix for L. The notion of the Weyl matrix is 
a generalization of the notion of the Weyl function (m-function) for the scalar case (see [T^ [2D]1 
and the notion of the Weyl matrix for the matrix Sturm-Liouville operator (see mm)- 

We will prove the following properties of spectral characteristics of L. 

Theorem 1. The pencil L has a countable set of eigenvalues; 

(i) they can he numbered (counting with their multiplicities) to satisfy the following asymp¬ 
totics as \n\ —)■ oo; 

Pnq = n + U}q + O {n~^) , u G w != {±0, ±1, ±2,... g=l,m, (:nge[0,l); (3) 

(a) for some sufficiently large N, all the eigenvalues pnq, \n\ > N, are real. 

The poles of the Weyl matrix M{p) coincide with the eigenvalues of L; for some sufficiently 
large N, 

(Hi) the poles at pnq, \n\ > N, are simple; 

(iv) the weight matrices anq '■= ReSp=p^^ ^(p) satisfy the asymptotics 

y ^ =—(I _hi)-^W^LW(I+ hi)-^+ 0(n-^) , Inl^cx), g = 1^, (4) 

mns Trn ^ ' 


and rankanq = ronq. Here mnq are the multiplicities of the eigenvalues pnq, Jq = {s-. Uq = Ws} 
and W,Iq are some m x m matrices (see notation before Lemma\^for details). 

Some of these properties were proved in [21] but for the matrix differential pencil with other 
boundary conditions, independent of the spectral parameter. Moreover, we need here more 
precise asymptotics for the recovery of the pencil. 

The properties, established in Theorem [T], are analogous to the ones for self-adjoint matrix 
Sturm-Liouville operators, but only for sufficiently large eigenvalues. The pencil can have a 
hnite number of nonreal eigenvalues and multiple poles of the Weyl matrix, even when the 
coefficients satisfy (II). This is possible even in the scalar case (m = 1). 

Example. Let m = 1, hi = LTi = ho = iLo = 0, Qi(x) = —fp, Qo(x) = q, where p and q 
are real constants. Clearly, conditions (I) and (II) are fulfilled. So we have a scalar boundary 
value problem 

y" + (P^ + 2pp + q)y = 0, y'{0) = yfir) = 0. 

Its eigenvalues are the roots of the equations p^ -|- 2pp + q = n^, n = 0,1, 2,.... Obviously, 
these equations can have nonreal roots. Now suppose p = g = 0. Then one can easily calculate 

COS pTT 

the Weyl function M(p) = -. Clearly, the eigenvalue p = 0 is a double pole of Mip). 

p sm pn 

For convenience, renumerate the eigenvalues {pnq}\n\<N by {Pn}n=TK- Denote by the 
multiplicity of pn as a pole of M(p) (then we assume pn = Pn+i = ••• = pn+m^-i). Let 
S' = {n G [1, iL]: n = lorp„_i ^ pn}. It is easy to show that the Weyl matrix admits the 
following representation: 


rrin — l 


M{p) = E E 


^n-\-u 


n&S u=0 


Pr. 


'li'+l 


(p 

|n|>Ar q=l 


^nq ^nq 


Pnq) 


( 5 ) 


where n = 1, iF, are some m x m matrices. □ 

^In case of multiple eigenvalues, the same weight matrices a„q will be counted multiple times. To avoid this 
problem, we divide the weight matrices in (|T]) and (O by the multiplicities rrinq. 
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We call the data SD := {p„, an}n=XK {Pnq, C(nq}\n\>N,q=i;rn spcctral data of the pencil 
L. This notion generalizes the notion of spectral data for scalar pencils (see m)- According to 
(15]) , the spectral data determine the Weyl matrix uniquely. Therefore, the inverse problem by 
the Weyl matrix, studied in im, is equivalent to the following one. 

Inverse Problem 1. Given spectral data SD, construct the coefficients of the pencil L. 

Thus, in this paper we develop a constructive solution of Inverse Problem 1. The analogous 
results for the scalar case were obtained in for pencils with arbitrary complex-valued 

potentials. However, in our paper we put an additional restriction (II) on the coefficients of the 
pencil, in order to guarantee the simplicity of the poles of the Weyl matrix at large eigenvalues 
and to work with asymptotics of the corresponding residues. The scalar case is easier from 
this point of view, because in the scalar case asymptotic formula (J3]) yields the separation of 
the eiganvalues for large n, so the poles of the Weyl function are automatically simple. In 
the matrix case, there can be inhnitely many groups of multiple eigenvalues, and the general 
situation is very difficult for investigation. 

The paper is organized as follows. In Section 2, we study the properties of the spectral 
characteristics of L, crucial for our method of recovery, and prove Theorem [H In Section 3, 
following the ideas of the method of spectral mappings (see [T91I22] ! we choose a model pencil L 
and then by integration in the complex plane of the spectral parameter, reduce Inverse Problem 
1 to a linear equation (see eq. (IT7|) and Lemma [ 8 ]). In Section 4, this equation is transformed 
into the so-called main equation (I57|) in a specially constructed Banach space. We investigate 
the solvability of this equation (Theorems [2] and [3]). As opposed to the similar equation for 
the Sturm-Liouville operator, the main equation for the pencil is uniquely solvable not for all 
X G [0, vr]. This difficulty appears even in the scalar case (see [7][8] ). In Section 5, we provide an 
algorithm for the solution of Inverse Problem 1 “step by step”, choosing a new model pencil at 
each step. We prove that our algorithm finds coefficients of L on the whole interval [0, tt] after a 
hnite number of steps. The algorithm, based on analogous ideas, was developed earlier for the 
matrix pencil on the half-line [16]. For simplicity, we work under the assumption A^ = 0, i.e. all 
the eigenvalues and weight matrices satisfy properties (ii), (iii) of Theorem [T] In Section 6, we 
show how to work with multiple poles of the Weyl matrix. We provide necessary modifications 
of the main equation and, finally, arrive at the solution of the inverse problem in the general 
case. 

Notation. Below along with L we consider pencils L = L{i, U, V) and L* = U*, V*): 

iY := Y" + {p^I + 2ipQi{x) + Qo{x))Y = 0, a; G (0, tt), 

U{Y) := r(0) + {iphi + ho)Y{t)) = 0, 

ViY) := Y’iir) + (ipH, + H,)Y{-k) = 0; 

Z := Z" Y Zt^p^I Y 2ipQi(px) Y Qo(x)) =0, x G (0, vr), 

U*(Z) := Z'(0) Y Z(0)(zphi Y ho) = 0, 

V*{Z) := Z'{7r) Y Z{TT){tpH^ + Hq) = 0, 

where Z is a row vector, and L*, formed in an obvious way. We agree that if a symbol 7 denotes 
an object related to L then 7 , 7 * and 7 * denote the corresponding objects related to L, L*, L*, 
respectively. We mention once more that the conjugate transpose is denoted by f. If Y (x, p) is 
a solution of ([ID, then Y*{x,p) := Y^{x,p) is a solution of tY* = 0. 

We consider the space of complex m-vectors C™' with the norm 

||y|| = max \yj\, Y = [yj]j=TYi^ 
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the space of complex m x m matrices with the corresponding indnced norm 

m 

11 11 ^ I I ’ ^ = Mj,k=T;^- 

k=i 

and the Hilbert space L 2 (( 0 , tt), C™) of m-vectors with entries from ^ 2 ( 0 , vr) with the following 
scalar product and the corresponding norm: 

Y\x)Z{x)dx= / yj{x)zj{x) dx. 

do 2=1 

{Y, Z) denotes the matrix Wronskian Y'Z — YZ'. 

In estimates and asymptotics, we use the same symbol C for different constants independent 
of X, p, etc, and r := Imp. 

2. Properties of spectral characteristics 

Let ip{x,p) = p)]j,fc=T;^ and S{x,p) = [Sjk{x, p)]j i.^j^ be the matrix solutions of 

equation ([T]) under the initial conditions 

s (0, p) = u (p) = 0, ^'(0, p) = p(x, p) = I. 

Introduce the matrix functions P+{x) and P_(x) as the solutions of the Cauchy problems 

F;(x) = ±Qi(x)P±(x), P±(0) = P ( 6 ) 

When condition (II) holds, they satisfy the relations 

P|(x)P+(x) = pI(x)P_(x) = /, X e [0,7r]. (7) 

It was proved in mi. that for X G [0, tt], z/ = 0,1, as |p| —)■ cxd: 

S^'^\x,p) = —exp(ipx)P_(x) + -—^ -exp(—ipx)P+(x) + 0(p^“^ exp(|r|x)), (8) 

ip^’'\x,p) = ■^^^exp(ipx)P_(x)(/ - hi) + exp(-ipx)P+(x)(/ + hi) + 

+ 0(p"-iexp(|r|x)). (9) 

Lemma 1. The pencil L has a countable set of eigenvalues. They can be numbered (counting 
with their multiplicities) in such a way that the eigenvalues have the following asymptotics as 

|?7,| —)■ OO.- 

pnq = n + Uq +O (n~^) , u G tc := {±0, ± 1 , ±2,... q = l,m, a;gG[0,l). 

Proof. It is clear that the eigenvalues of the pencil L coincide with the zeros of the characteristic 
function detH(p). This function is entire and has an at most countable set of zeros. By fl9l) . 

l/(p) = ip|^exp(ip 7 r)(J + Pi)P_( 7 r)(J - hi) 

- ^exp(-ip 7 r)(/-Pi)P+( 7 r)(J + hi)| + 0 (^exp(|r| 7 r)j, (10) 
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( . \ m 

j det -P-(vr) det (/ + Hi) det(exp(ip7r)J — exp(—ip7r)y4) 

+ 0^|pP“^exp(m|r|7r) |p|-)■ co, 

where 

A = Pz\ti) (/ + Hi)-^ (/ - Hi) P+(7r) (/ + hi) (J - hi)-'. (11) 

It follows from ([7]) and (II), that detP_(7r), det (/+ Pi) and det (J —Pi) are nonzero. 
Therefore the eigenvalnes coincide with the zeros of the analytic fnnction of the form 

A(p) ;= det^exp(2p7r)J — exp(—ipTr)^ + exp(|r|7r)) j, |p| —)■ cxd. (12) 

The fnnction Ao(p) := det(exp(ip7r)/ — exp(—zpTr)^) has zeros 

plq = n + ujq, n e Z, q = 

where ^ Inpg, Pg are the eigenvalnes of the matrix A. 

Since by ([7]) and (II) the matrices P-(7r), P+(7r), (/ + Pi)”' (/ — Pi) and (J + hi) (/ — hi)”' 
are nnitary, the matrix A is also nnitary. Its eigenvalues are simple and have the form pg = , 

where ipq can be chosen in the interval [0, 27r). Consequently, cUg = ^ G [0,1). 

Consider the circles 7 „g with the centers and radii equal to h, such that 

0 < h < ^ inf |p°g - p1^p\ , 

where the inhmum is taken over all n, h G Z, 1 < g,p < m, such that p^^ ^ 

Note that |Ao(p)| is bounded from below for p G 7„q, n E Z, q = l,m. Furthermore, 
|A(p) - Ao(p)| = 0{n~^) yields |A(p) - Ao(p)| < |Ao(p)|, for p G 'fnq, n G Z, g = l,m for 
sufficiently large |?7,|. Therefore we can apply Rouche’s theorem to the analytic functions Ao(p) 
and A(p) — Ao(p) in the circle 7 „g and conclude that A(p) has the same number of zeros as 
Ao(p) in this circle. As 6 tends to zero, we get 

Pnq= Pnq+^q + ^nq, £ng = o(l), |n|-)-CX), q = 1,171. (13) 

Since the matrix A is unitary, it can be represented in the form 

A = W^MW, ITt = Vp-i, M = diag{pi,p2,...,Pm}. (14) 

Substituting flT^ and flTTD into flTT]) . we get 

det ^p.g exp(27rP„g)/ — Ai + O (’^-') j = 0. 

Since Ai is diagonal, we conclude that 

£nq = O (n,-') , \n\ —)■ oo. (15) 

Consider the contour T^^ = {p: |p| = R} with a sufficiently large R, that does not contain 
points p°q. Apply Rouche’s theorem to it to the functions (ip)'"A(p) and (ip)”’'Ao(p). As a 
result, detR((p) has the same number of zeros as (ip)™'Ao(p) inside this contour, and together 
with 0131) and 0151) this yields the assertion of the Lemma. □ 
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Lemma 2. The following relation holds 


j ip^x,Png)(p{x,Pnq)dx = ^ (^I + h\hi'^ + O {n ^) , |n|-)■ CX), q = l,m. (16) 

0 

Note that (^I + h\hi^ >0. 

Proof. We will use the notation [/]=/ + 0{\p\~^). According to ([9]), 

p\x,p)p{x,p) = ^|^exp(-ipa;) (^I - h\^ P[{x)[I] + exp{ipx) (^/+ hjj P|(a;)[J]^ 

■ {exp{ipx)P_{x) (/ - hi) [I] + exp{-ipx)P+{x) (/ + hi) [/]). 

Using ([7]), we derive 

p\x, p)p{x, p) = ^(^I + hlhi^ [I] + ^ exp(2zpa;) (J - h) pI{x)P_{x) (J - hi) [I] 

+ ^exp(-2ipa;) (/ + hi) P}_{x)P+{x) (/ + hi) [/]. (17) 

Integrating by parts, we obtain 


j exp{2ip7i)Pl{x)P.{x) dx 
0 


1 

2ip 


exp{2ipx)Pl{x)P-{x) 


TT 

0 


1 

2ip 


J exp {2ipx)Pl{x)QI {x)P-{x) dx 
0 


0{n ^) , 


P Pnqi 


and analogously, 

J expi-2tpx)Plix)P+{x) dx = 0 (n-i) , p = p^g. 

0 

Now it is clear that flTBD follows from f[T7)l . □ 

Lemma 3. Starting from some number N , all the eigenvalues png, \n\ > N, q = l,m, are real. 

Proof. Let po he an eigenvalue and Y{x, po) be a corresponding eigenvector for the pencil L. 
Consider the scalar product in L 2 (( 0 , vr), C™): 

5[F] := (y, Y" + [pII + 2*poQi + Qo) = {Y, Y") + pl{Y, Y) + 2tpo{Y, QiY) + {Y, QoY) = 0. 
Integrating by parts, we get 

(F, Y") = -Y^{7r){tpoHi + Ho)Y{7r) + Y\0){tpohi + ho)Y{0) - (W, Y'). 


Every eigenvector can be represented in the form 


Y {x, Po) = p{x, po)Cg,{po), ( 18 ) 

where C<^(po) is a constant column vector, it can be chosen so that ||C^(po)|| = 1- 
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It follows from ([9]) and ffTSjl . that when |po| —t cxd by the set of eigenvalnes, 


Y(x,p„) = 0(1), 


( 19 ) 


nniformly with respect to x G [0, tt]. 
Let po = a + ir, r ^ 0. Then 


Im 5[F] = 2aT(V, V) + 2t(V, QiV) - tY^{T r)iHiY{Tr) + tY^{ 0)ihiY{0) = 0, 


a = 


2(Y, QiY) - Y\7i)iHiY{n) + F^(0)2hiF(0) 


( 20 ) 


It follows from ffT^ . that the nnmerator of this fraction is bonnded. Lemma [2] helps us to prove 
that the denominator is bounded from below. Indeed, by positiveness of the matrix {I + h\hi), 
we have 


{X,{I+ h\hi)X)>a\\Xf, a>0. 


Hence 

TT 

{y^y) = j Clp\x,po)p{x,p)C^dx > a\\C^f 

0 

for sufficiently large |po|- 

In fact, we have proved that the real parts a of eigenvalues are bounded as |po| —t oo. But 
this contradicts the asymptics from Lemma [TJ Thus, r = 0 for all eigenvalues, sufficiently large 
by their absolute value. □ 


Remark 1. Suppose 


F[Y] := iY\ Y') - (F, Q^Y) + Y\n)H^Y{Ti) - F+(0)hoF(0) > 0 


for all eigenvectors Y = F(x, po)- For example, this holds when Qo(x) < 0, iLo > 0, ho < 0. 
Let G\Y] be the numerator of the fraction in fl20|) . Then 

ReS[Y] = {a^ - t^){Y,Y) + aG[Y] - F[Y] = 0. 

Note that by (12(11) . crG[F] = —2a‘^{Y,Y). Hence 

F[Y] = -(a^ + t^)(Y,Y) < 0. 

We conclude = 0, so in this special case, all eigenvalues are real. Moreover, we note 

that fl20|) can hold only for po = 0. 

It is easy to check that 


$(x,p) = ^(x,p) + (p(x,p)M(p), 

(21) 

M(p) = -(F()p))-V(^). 

(22) 


Consequently, M(p) is a meromorphic function and its poles coincide with the zeros of det V (<p), 
i.e. the eigenvalues {pnq}- 


Lemma 4. The poles of the matrix function M{p) at p = pnq for sufficiently large |n| are 
simple. 
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Proof. We will prove that the matrix function {y{(p)) ^ has only simple poles at sufficiently 
large eigenvalues, using the following well-known fact (see m Lemma 2.2.1]): 

If there do not exist two nonzero vectors a and h, such that 


V{(p)a = 0, 


(23) 


at some point po, then po is a simple pole of {y{(p))~^ 

In order to use this fact, we derive a formula for 

In this proof, we assume that p is real. Let 'ip{x,p) = [ipjkix, p)]j^j^^i^ be a matrix solution 
of equation ([T]) under the initial conditions 


V’(7r,p) = /, I/(V>) = 0. 

In view of (II), fj^x^p) is a solution of the following equation 

Z" + Z{p^I + 2ipQ^{x) + Qo{x)) = 0, (24) 

and consequently, the Wronskian p),Lp{x, p)) does not depend on x. Therefore we calcu¬ 

late 

{'il)\x,p),p{x,p)) = {'il)\x,p),p{x,p))\^=o = {U{'ll)))\ 

{'ip\x,p),p{x,p)) = {f)\x,p),Lp{x,p))\^=^ = -Id(<p), 

= -V{p). (25) 

Let Po be a real eigenvalue of L. Then 

^(V’^(x,p),(p(a:,po)) = {po - p)'iIj\x, p){{p + po)I + 2iQi{x))ip{x, Po), 

{'il)\x, p), (p(x, Po)) Ip = (po - p)Lfi(p(vr, Po) 

- I/(p(x,po)) - (u^ifix^p))^ - (po -p)V’H0,p)hi. 

By fl25|l we have 


^ p)) - V (p(x, Po)) = (po - p) 


V''^(0,p)hi - iLi(p(7r,po) 


-f- / 'ip^x, p)((p + po)I + 2iQi(x))ip(x, Po) dx 


Taking a limit as p tends to po, we get 


A 

dp 


^{p)\p=PO 


V'^(0,po)hi - iLip(7r,po) 


+ 2 


J ^jJ\x,po){poI+ iQi{x))p{x,po)dx 


(26) 


0 

Every eigenvector Y{x,po), corresponding to the eigenvalue po, has the form 

F(x,po) = p(x,po)C<^ = fj{x,po)C^, 

where and are some constant vectors. We can assume that ||C^|| = 1. Using (126]) . we 
obtain 






( 27 ) 




TT TT 

+ 2 / y^(a:,p)iQiF(a;,p)rfa; + 2po / C'J,(^'^(a:, po)‘^(a;, 


If ||C'<^|| = 1, then Y{x,po) = 0(1), |po| ^ C )0 and Cp = (p(7r,po)C'^ is also bounded. By 
Lemma [21 the integral in the last term is positive and bounded from below for sufficiently large 
IPoI- Therefore, C]p^V{p)C^ tends to ±cxd as po = Pnqi n —)■ ±oo. 

If a and b satisfy fi22]l . then a gives an eigenvector, i.e. a = C^. On one hand, from the 
second relation of fl23|l we get 

clv(p)b = -cl-^v(p)c,. 

that is unbounded for sufficiently large |po|. On the other hand, by (I2BD we have C^V(p) = 

—O^(t/('0))^= 0. We arrive at a contradiction, so for sufficiently large |po|, the vectors a and 
b in fl25D must be zero and the poles of {V{ip))~^ are simple. □ 

Remark 2. Consider the case from Remark [H Note that if po is real, then the right hand side 
of (177)) equals zero if and only if (I2UD holds. But according to Remark [H fl2UD can be valid only 
for Po = 0. Moreover, if F[y] > 0 for all eigenvectors belonging to the zero eigenvalue, then 
(1201) never holds. Following the proof of Lemma IH we conclude that in this special case all the 
poles of M(p) are simple (if F\Y] can be zero for po = 0, then all the poles except po = 0 are 
simple). 

Introduce the notation: Jq = {s\ ojs = Wg}, Iq = 


— 



j = k and ooj E Jq] 
otherwise. 


Let uinq be the multiplicity of the eigenvalue pnq (as a zero of the analytic function det V ((p)). 
We consider the data anq ■= Res M(p) and m^q only for sufficiently large |n|, for which the 

P=p7iq 

assertions of Lemmas ITlid] and 01 hold, so here we do not deal with non-simple poles. 

Lemma 5. The ranks of the weight matrices anq are equal to mnq for all such n, q that the 
corresponding poles of M{p) are simple. 

The proof is similar to the one for the matrix Sturm-Liouville operator [T01 Lemma 4]. 

In the following lemma, we split the eigenvalues into groups according to their asymptotics 
and sum all residues corresponding to a group. 

Lemma 6. The following relation holds 

y ^ =—(I-hp-^W^LW(I+ hi)-^+ 0(n-^) , Inl^cx), g = 

mns Trn ^ ^ 

S^Jq 

where W is the unitary matrix, defined in (HI, 

Proof. Using (jS]), we obtain 

y{S) = ^exp(2p7r)(J + Hi)P_{7r) + ^ exp(-ip7r)(/ - i/i)P+(7r) + O (^|p|"^ exp(|r|7r)), 
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Substitute this and flTOj) into fl2^ : 


M(p) = -(l/(v.))-V(S) = -M-‘|lexp(v^)(/ + H,)P_(r)(I - h,) 

- 2 - ^^i)^+(7r)(/ + hi) + o(^|p|"^exp(|r|7r))| 

• exp(ip7r)(/ + Hi)P_{n) + ^ exp(-zp7r)(J - Hi)P+{7i) 

+ 0(^|p|"^exp(|r|7r) = -{ip)~^{I - hi)"^|exp(2ip7r)/ 

- Pl\7r){I + - i/i)P+(7r)(/ + h){I - h)-^ + 0(|p|-' exp(2|r|7r)) I 

• |exp(2ip7r)/ + Pr^(vr)(J + - Pi)P+(7r) + 0(^|p|"^ exp(2|r|7r) j |. 

Finally, we get 

M(p) = - hi)“^|exp(2ip7r)/ - A + O(^|p|“^ exp(2|r|7r) j | 

■ |exp(2ip7r)J + B + 0^|p|“^ exp(2|r|7r) j (28) 

where A was dehned by flTTD and 

B = A{I - hi){I + h)-\ (29) 


Introduce the matrix function 

Mq{p) = —{ip)~^{I — hi)“^(exp(2ip7r)/ — A^ ^ ■ (exp(2ip7r)/ + P). 

Let 'ynq be the contours used in the proof of Lemma 1. They are circles with centers in 
Png = n + ujq and a sufficiently small fixed radius, such that p'^g G int 'jnq if and only if 

. 

rms rnq 

According to (l28|) . we have 

||M(p) - Mo(p)|| < ^, pE'Jnq- 

The residue theorem yields 

-—: [ (M(p) — Mq(p)) dp = - Res M(p) — Res Mo(p) = 0(n~‘^), \n\ —)■ oo. (30) 

2:711 J T^ns p—pnq P—Pnq 

luq 

Clearly, Mq{p) has simple poles at p^^, and using flTTD . we derive 

Res Mo(p) = — lim (p — po)(ip)“^(/— hi)“^IF^(exp(2fp7r)/— Ad) ^IF(exp(2fp7r)/+ P) 
p=p% p^p% 

= - h)-'W'ppI,W{p,I + B) =: al,. 

^^Pnq 

Using the obvious relation p~^IqA4 = Iq and ([29]), we get 

W^p-^IqWB{I + hi) = W^p-^IqWW^MW{I - hi) = WUqW{I - hi). 
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Therefore 


1 




■{W^IgW{l + hi) + W^IgW{l - hi)), 


nq 


a. 


nq 


= —{I - hi)-^fT%fT(J + hi)-i + 0(n-2), 


\n\ —)■ oo. 


nn 


Together with fl5U]l this yields the assertion of the Lemma. 


□ 


Summarizing the results of the Lemmas, we arrive at Theorem [H 

Note that there is a sufficiently wide class of pencils, for which assertions (ii) and (iii) of 
Theorem [T] are valid for all n E w. Namely, (ii) holds for all eigenvalues when 


F[Y] := (r, r) - (F, QoY) + Y\7r)HoY{n) - Y^{0)hoY{0) > 0 


for all eigenvectors Y = Y(x,po). For example, this holds when Qo(x) < 0, ffg > 0, ho < 0. 
In this case, (iii) is valid for all eigenvalues except zero. To make it valid for zero, we require 
additionally that F\Y] > 0 for eigenvectors corresponding to po = 0. 

3. Contour integration 

Assume that the spectral data SD of an unknown pencil L are given. 

At first, choose a model pencil L with the spectral data satisfying asymptotics ([3]) and (jl]) 
with the same constants as the data of L. This can be done by the following algorithm. 

Algorithm 1. Choice of a model pencil. 


1. Construct M(p) via (jS]). Using (12511 . one can easily derive 

M{p) = {ip)~^{I + hi)~^ + p = ir, r ^+oo. 

Therefore, having M{p), we can obtain hi and set hi: = hi. 

2. Find ojq from the asymptotics ([3]), Pq\ = M. := diag{/iq}g^i^. 

3. Construct Jq and Iq. From the asymptotics (0]), using hi, find W^IqW . Note that 

m 

A = W^MW = \Jq\-^PqWUqW, 

9=1 


SO we can construct A = A. 

4. Put Hi = 0. According to flTTll . we need to find such a potential Qi{x), that 

Pl(7r)P+(7r) = A{I - hi)(/ + hi)-i =: C. 

Try a constant anti-Hermitian matrix: 

Qi(a;)=T, pt =-T. 


Then P±{x) = exp(±Ta;), and we arrive at the matrix equation exp(2T7r) = O. It can 
be easily solved, since the matrix O is unitary. 

5. Put Qoix) = ho = Hq = 0, X E [0, tt]. 
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Denote 


(31) 


D{x,p,e) = 


{p*{x,e),p{x,p)) 

p-e 


Similarly to the scalar case (see [231 Lemma 3]), one can obtain the following estimate 

\\D{x,p,e)\\ < ^ exp(|Imp|a;)exp(|Img|a:), (32) 

\p-e\ + i 

where Cx is a constant depending on x. 

Consider the block-matrix V{x,p) = [Pjk{x, p)]j^k=i ,2 dehned by 


V(x, p) 


p{x,p) 

$(x,p) 


p{x,p) 

$(a;,p) 

_ (p'{x,p) 

^\x,p) _ 


_ t'{x,p) 

<h'(a;,p) _ 


(33) 


Note that if Y satishes ([T]) and Z satishes fl2^ for the same p, then 

L(z.r) = o. (34) 

Hence the Wronskian {Z, Y) does not depend on x. Conseqnently, nsing bonndary conditions 
for (/?, (p*, $ and $*, we easily derive 


p{x,p) 

$(x,p) 

-1 

1 

O 

* 

1 

* 

_ 1 

_ t'{x,p) 

$'(a:,p) _ 


. -T*'{x,p) (p*(x,p) 


Vji{x,p) = p^^ ^\x,p)¥'{x,p) ^\x,p)^*'{x,p), 

Vj2{x,p) = <h(^"^)(a;,p)(p*(a:,p) - p)^*{x, p). 

We mention two other facts that follows from 0341) and will be used further: 

{p*{x,p),p{x,p)) = {p*{x,p),p{x,p))x=o = 0, 
(4>*(x,p),$(x,p)) = M{p)-M*{p) = 0, ^ M(p) ^ M*{p). 

Denote 


(35) 

(36) 

(37) 

(38) 


Gs := {p e C: Ip - Pnq\ > ^, Ip - p^gl > 6, n e w, q = l,m}, (5 > 0. 

In a standard way (see [IS]), one can prove the following Lemma. 

Lemma 7. Assume hi = hi. Then for each fixed x G [0, tt], the entries of the block-matrix V 
have the following asymptotics as |p| oo, p G Gs-' 

Vii{x,p) = Q{x) + 0{p-^), Vi 2 {x,p) = p-^A{x) + 0(p"^), 

V 2 i{x,p) = -pA{x) + 0(1), V 22 {x,p) = D(x) -h 0(p"^), 


where 


kl{x)-.= - {P.{x)Pl{x) + P+{x)Pl{x )) , 

A(a;) 1 (p_(i)p!(i) - P+{x)pI(x) ) . 


(39) 

(40) 


For simplicity, we assume that = 0 and iV = 0, i.e. the assertion of Theorem [T] holds for 
all the eigenvalues of both pencils L and L. We return to the general case in Section 6, since 
it requires some modihcations. 


12 













Lemma 8. The following relations hold: 


m 

n{x)(p{x,p) = (p{x,p) + Png)m-gar,gD{x, p, Png) - (fix, Png)fh-gangD{x, p, Png)y 


n£w q=l 


(41) 


m 

D{x, p, 9) - D{x, p, 9) - p*{x, 9)A{x)(p{x, p) + ^ ^|L>(a;, pnq, 9)m~lanqD{x, p, pnq) 


n£w q=l 


D{X, Pnq, 9)m^^anqD{x, p, fnq) \ = 0. (42) 


Proof. 1. Choose a number 6 > supmax{|Imp„g|, |Imp„q|}. Denote 


n,q 


Eh := {p: |Imp| < 6, IpI < i?}, jr := OEr, Th := {p: |p| = i?}, 7 % = TR-'yR 

(the dehned contours have counterclockwise circuit). By Cauchy’s integral formula, 


Vik{x,p) - ^{x)6ik = — 


Vik{x,9) - Q{x)6ik 




P-9 


d9, /c = l,2, xG[0,7r], pGinty^, 


where Sjk is the Kronecker delta. Put 

£r{x,p) : = 


'Pik{x,9) - Q{x)6ik 
p-9 


d9, 


Since we have chosen hi = hi, we can apply Lemma [7] and obtain lim eR{x,p) = 0. Conse- 

R-^oo 

rR<zGs 

quently, 


Vik{x,p) = VL{x)6ik+ lim — 

R—^oo li'KX 
Vr^Gs 


V\k{x,9) 
p-9 


d9, A; = 1,2, xG[0,7r], |Imp| > 6. (43) 


Substituting (1451) into (1551) and using (1551) . (15T]) and (15T]) . we obtain 

1 


ip{x,p) = n{x)(f{x,p) - lim . 

ri-sa X, 


p(x, 9){M{9) - M{9))D{x, p, 9) d9, 


since the terms with S{x, 9) and S*{x, 9) vanish by Cauchy’s theorem. Calculating the integral 
over 7 ij by the residue theorem, we arrive at (TO- 

2. Analogously to (1451) . for x G [0,7r], |Imp|, |Im6*| > b and jr C Gs, we derive 


V2i{x,p) -P2 i{x,9) a / n , 1 

-= -A(a;) + 


r2i{x,o 


p-9 

Vjk{x,p)-Vjk{x,9) _ 1 


2vri (p-Oi^-^) 
'Pjk{x,f) 


df + eR^2i{x,p,9), 


p-9 


27ri (p-0(^-^) 


+ £Rjk{x, p, 9), (j, k) ^ (2,1), 


where hm eRjk{x, p, 9) = 0, j, k = 1, 2. 

R-^oa 

TrCGs 

For an arbitrary vector function Y = Y{x), in view of (13 5 p and (133 p . we have 


r 


Y 

Y' 


T 

P' 




<F 


{(p\Y). 
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Therefore 




V{x,p)-P{x,9) 

p-e 


p{x,p) 

(p\x,p) 


= p>*{x,e)A{x)(p{x,p) 


1 f {{^*{x,d),^{x,0){^*{x,^), p{x,p)) - {p*{x,e),^{x,^)){p*{x,^),^{x,p)))d^ 
27ii L, 


{p-0(^-0) 


+ e]i{x,p,e). (44) 


Using fl3^ and fl35|l . we obtain 

[p*'{x,e),-p*{x,e)]'P{x,p) 


p{x,p) 

p'{x,p) 


[p*'{x,e),-p*{x,e)]v{x,e) 


(p{x,p) 
(p'{x,p) 

Combining fl44)) . fITS]) and fITBD and using fl^ . we get 


= {p*{x,e),p{x,p)), 
= {(p*{x,e),(p{x,p)). 


(45) 

(46) 


D{x,p,e) - D{x,p,e) 


(p*{x,9)A{x)p{x,p) 


-^/ [d{x,^,9)M*{0D{x,p,0 

- D{x, 9)M{^)D{x, p, O) d). 


Note that M(p) = M*[p) by fl55]l and limjj_>.o £\j{x, p, 9) = 0. Therefore, calculating the integral 
over 7 jj by the residue theorem, we obtain (02]). □ 

Denote 

D(a;) ;= ^ [P-{x)P^_{x) + P+(a;)P|(a;)j , 

A(a;) := ^ (p_(a;)Pl(a:) - P+(x)Pi(a:)) , 

(compare with fl39l) . fHOjl h Symmetrically with the relations flTT|) and f02|) . we have 


m 

D(x)v7(x,p) = (p{x,p) -^^^(p{x, Pnq)m~lanqD{x, p, Pnq) - (p{x, Pnq)m~lanqD{x, p, pnq)^ , 

n£w q=l 

(47) 


D(x,p,e) - D(x,p,9) + <p’(x,S)A(x)ifi(x,p) + ^ ^|D(a:,p„„9)m„y„,D(a:,p,p„,) 

n£w q=l 

- D{X, Pnq,9)ThPganqD{x, p, Pnq)^ = 0. (48) 

4. Construction of the main equation in a Banach space 

In the previous section. Inverse Problem 1 was reduced to the linear equation flTT|l . However, 
the series in this equation does not converge absolutely, so it is inconvenient to use flTT|l as a 
main equation of the inverse problem. Now we plan to transform flTTl) into an equation in a 
Banach space of infinite sequences. 

In contrast to the scalar case, the spectrum of the matrix pencil L can contain subsequences 
of eigenvalues with the similar asymptotics (some of cOq in (|3D can be equal). In order to 
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overcome this difficulty, we divide the eigenvalues into groups Gk by asymptotics and then 
operate not with single eigenvalues but with the whole groups. 

In this section, we assume that x G [0, tt] is fixed. 

Introduce sets Gnq = {Pns,Pns, s G Jgj, n G w. If some ujq are equal to each other, the 
corresponding sets coincide. Therefore we assume (without loss of generality), that 
are all the different values among Consider Gnq, q = I,p and renumerate them as 

Gk, k & Z\{0}, in the natural order. 

Let G = {giYi=i be a non-empty set of distinct numbers. Define the diameter of G\ 

diam(G) := max {\gi - fi-i |} 

z=2,r 

and a finite-dimensional space B{G) := {/: G —>■ with the norm 

||/||b(g) = maxj |l/(^i)||, maxj ||/(^i) - f{gi) 

Note that diam(G) can be zero only in the case r = 1, but then ||/||b(g) = ll/(5'i)ll- 
Consider a Banach space 

B = {f = {fk}kez\{o}- fk ^ B{Gk), 11/11^:= sup ||/fc||B(GA,) < co}, (49) 

kez\{o} 



By (13]), diam(G„„) < —r, and consequently, diam(Gfc) < 

|n| 

Lemma, it is easy to see that 


W\' 


Using dnj) and Schwarz’s 


\\p{x,pi)\\ < G, \\p{x,pi) -p{x,p 2 )\\ < G|pi -P 2 I, 
for pi,p 2 G Gfc. Therefore ip{x,p) forms an element of B: 

ip(x,p)lB = WY,P)\GnAneu.,q=T:^ ^ ^ 

P(x^P)\Gn, = W(x,Pns), (p{x,Pns), S G Jq}. 

Denote ■^(x) := (p(x,p)|s, YY) •= vYyP)\B- Then the relations flTT|l and fl47l) can be trans¬ 
formed into the following equations in the Banach space B: 

VL{x)'ii){x) = ip{x){I + R{x)), (50) 

D(x)V’(a:) = '?/’(x)(/ — R{x)), 

where I is the identity operator in B, and R{x), R{x) are linear operators, acting from B to B. 
[§ The explicit form of R{x), R{x) can be derived from flTTl) and flT7|) . Further we investigate 
the operator R{x), the same properties for i?(x) can be obtained symmetrically. 

According to (jUj), the operator R{x) acts to an arbitrary element ip ^ B in the following 
way: 

(V’i?(x))^ = ^ ipkRk,nY), Rk,nY)- ^{Gk) B{Gn), k,n & Z\{0}. (51) 

| fc |>0 

Lemma 9. The series in fl5T|) converge absolutely and the operator R{x) is bounded on B. 

^These operators belong to the ring for which the elements of B form a right module, so we write operators 
to the right of operands. 
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Proof. Consider two sets G„ = {giYi=i and Gk = {sj}Yi The operator P) in B{Gn) 

in the following way: 


Hai) = fRk,n(x) = 'P^f(aj)ajD(x,gt,Sj). 


(52) 


i=i 


where Uj = mYocnq for Qj = pnq and aj = -m^^anq for gj = pnqM 
Let us prove that 


heBiGn), \\h\\BiG„)< 


C{\n\ + \k\) 
k'^{\n — /c| + 1) 


ll/l|(!(Gi), 


(53) 


Using fl^ . we get 

^(^i) = /(si) ) D{x,gi,si) + f{si)'^aj{D{x,gi,Sj) - D{x,gi,Si)) 

i 

+ “ fYi))ajD{x,guSj), 


0=1 


i=2 


h{gi) - h{gi) = Y1 9i, Sj) - D{x, gi, sj)) 


i=2 


i=i 


/(^i) iD{x, gi, si) - D{x, gi, si)) + /(si) ^ {D{x, gi, Sj) - D{x, gi, Sj) 


0 = 1 


i=2 


- D{x, gi, si) + D{x, gi, si)) + E(/( Sj) - /(si))«i {D{x, gi, Sj) - D{x, gi, sj)). 

i=2 

Since / G B{Gk), we have 

||/(s.)l| < ll/ll8(G.). ||/(i>j) - /(i>i)|| < |A^|-‘|I/I|b(g.). 

Using dS]), dl]), fl5^ and Schwarz’s lemma (see USD. it is easy to obtain the following estimates: 

i 


a. 


G „ „ G 

- p’ - jky 


\\D{x,gi,Sj)\\ < 


E 

i=i 

C{\n\ + \k\) 
(|n - fc| + 1)’ 


\D{x,gi,Sj) - D{x,gi,Sj)\\ < 


\D{x,gi,Sj) - D{x,gi,si)\\ < 
G{\n\ + \k\) 


C{\n\ + \k\) 
\k\{\n - fc| + 1)’ 


|n|(|n — k\ + 1)' 


\D{x,gi,Sj) - D{x,gi,Sj) - D{x,gi,si) + D{x,gi,si)\\ < 


G{\n\ + \k\) 


|n||fc|(|n — fc| + 1) ’ 

Combining these relations and using the dehnition of S(G„)^rLorm, we arrive at fl53|l . 
Substitute flS5|) into fl5T]) : 


||'0-R(a;)|| = sup {YR{x)) 


nGlj 


< 


E 

Jfc|>0 


\n\ + \k\ 


k‘^{\n — k\ + 1) 


<CM- 


Hence ||-R(a;)| 




< oo. 


□ 


^Here, in order to represent the operator in a more convenient way, we consider G„ as a multiset. For 
example, if Pnq = Pnq, they correspond to two equal elements in G„: gi and gj, f{gi) = figj). Clearly, this 


does not affect the norm ||/||e(G„). 
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Note that fHHjl is equivalent to the relation 

R{x) — R{x) + 0(x) + R{x)R{x) = 0, (54) 

where the operator 0(a;): B ^ B for each hxed x G [0, vr] is dehned by the following formula: 


ipQ{x) = 


nq )m^lanq^*{x,Pnq) “ 'lp{pnq)m^lanq<^* {x, Pnq)] 

n£w q=l 


A{x)p{x,p)\B. (55) 


Theorem 2. For each fixed x G [0,7r], the subspace of solutions of the homogeneous equation 

/3(J + ^(x)) = 0, /3 gS, (56) 

has dimension m — rankQ{x). In the particular case detf2(x) 7 ^ 0, equation fl56l) has only the 
trivial solution. 


Proof Let (3 G B he a solution of fl56|l . Then by fl541i . 

fi{I + R{x)){I - R{x)) =13- I3e{x) = f3- K{x)ij{x) = 0, 

where 'fi{x) = p{x,p)\B and K{x) is a row vector, which can be found from fl5^ . Using fl5U]l . 
we get 

I3{I + R{x)) = K{x)fi{x){I + R{x)) = K{x)n{x)fi{x) = 0. 

Let IC{x) be the space of solutions of the equation K{x)Q{x) = 0 for each hxed x G [0, tt]. It is 
clear that dim/C(a:) = m — rankf2(a:), and /3 is a solution of fl56|l if and only ii (3 = K{x)fi{x), 
K{x) G IC{x), this yields the assertion of the Theorem. □ 

Theorem 3. For each fixed x G [ 0 , 7 r], such that detff(a;) 7 ^ 0, equation 

z{x){I + R{x)) = •ip{x), (57) 

has a unique solution z(x) = (Q(x))~^if(x) in the Banach space B. 

Equation fl571) is called the main equation of Inverse Problem 1. 

Remark 3. For the matrix Sturm-Liouville operator (when Qfix) = 0), we have f2(x) = /, 
A(a;) = 0. Therefore, the main equation is uniquely solvable for all x in [0, tt] and by flHTjl the 
operators (/ + R{x)) and (/ + R{x)) are inverses of each other. 

Remark 4. The operators R{x) and R{x) can be approximated by sequences of hnite- 
dimensional operators in B. Indeed, let Rl^ = Rk,n for all n G Z, 1 < \k\ < s, and all 
the other components of R^ be equal to zero. Using the estimates for Rk^n from the proof of 
Lemma [9l it is easy the show that lim \\R^ — R\\b^b = 0. 

S—>-oo 

Therefore the operator R{x) is compact, and one can apply the hrst Fredholm theorem to 
the operator (/ + R{x)) and obtain from Theorem [2] even more general result than Theorem [3l 
if det fl{x) 7 ^ 0, the operator (/ + R{x)) has a bounded inverse on B. 


5. Algorithm for the solution 

In this section, using the main equation fl57jl and Theorem [3l we provide a constructive 
algorithm for the solution of Inverse Problem 1. 

Let X G [0, vr] be a hxed point, such that det fl{x) 7 ^ 0. Suppose we have solved equation 
m and found z{x). Multiplying (ITT]) by (Tl(a:)) we get a formula to construct z{x,p) : = 
{Q{x))~^(p{x, p) from z{x) = z{x,p)\B- 


17 




In a similar way to Lemma [HI it can be proved that 


f2(a:)<I)(a;, p) = <h(a:,p) 


m 


n£w q=l 


^nq ^nq 


E{X, p, Pnq) - Pnq)m^l ^nqE(^X, p, 


E{x,p,9) : = 


{p*{x,e),^{x,p)) 

p-e 


Therefore, knowing z{x, p), we can build the matrix function w{x, p) := (hl(a:)) ^<l>(x, p) by the 
formula 

m 

W[X,p) = <h(a;,p) Png)m-ganqE{x, p, Pnq) - z{x, Pnq)fh-gangE{x, p, Png)y (58) 


n£w q=l 


Denote z*{x,p) := z^{x,p), w*{x,p) := w^{x,p). Then ip*{x,p) = p'^(x,p) = p)!!^(x), 

$*(x, p) = <l)'^(a:,p) = w*{x, p)fl^(x). Using flHH]) . we get 

0 = (p$* — <l>(p* = ^l{zw* — 

I = (p<l)*' — ^{p*' = Q[zw* — wz*)Q^' + Q{zw*' — wz*')Q^. 

Therefore, 

zw* - wz* = 0, 

{zw*' — wz*')~^ = (59) 

It follows from flH7)l that (2;*Dl,r2^) = 0. Consequently, 

= z*Q^Qz' - z*'Q^Qz. (60) 

If the functions z, z*, w and w* are known, one can hnd from fl59l) . Then the right- 
hand side of fionj) is known, and we can calculate Dl D — DlD'. Differentiating we hnd 
D -|- DlD'. Then we can hnd and calculate 


a{x) := = (D^D)-^DtD'. 


Solving the Cauchy problem 

D'(a;) = a(a;)D(a;), D(0) = J, 

we get D(x) and then can hnd p{x,p). 

We arrive at the following algorithm for the solution of Inverse Problem 1. 

Algorithm 2. Let the spectral data of the pencil L be given. 

1. Choose a model pencil L, using Algorithm 1. 

2. Construct the space B, associated with the spectral data of L and L, 'ip{x) and R{x). 

3. Find z{x), solving the main equation fl57)l . 

4. Construct z{x,p) and w{x,p), using fjTTjl and fl58l) . then get Dl(x)D(x) via fl59ll . 

5. Find D(a;) and p{x,p) = Q{x)z{x, p). 

6. Substituting p{x,p) into ([I]) and ([2]), obtain Qs{x), hg, Hg, s = 1,2. 
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Note that Algorithm 2 works only for the case detf2(x) 7 ^ 0 for all x G [ 0 , 7 r]. Further we 
provide a modihcation of this algorithm for the general case. 

First let us modify step 4 of Algorithm 1 to construct a pencil L with spectral data satisfying 
([3]) and dl]) (with the same constants as the spectral data of L) and moreover, satisfying the 
condition 

Qi{x) = Qi{x), 0 < X < 6 < TV. (61) 

Note that the equation 

iI + H,)-\l-H^) = 0, (62) 

where the matrix O is unitary, has a solution 

H, = {I + 0)-\l-0), (63) 

if and only if det(/+C>) 7 ^ 0. In this case, the matrix Hi satishes (I) and (II), i.e. det{I±Hi) 7 ^ 0 
and Hi = -HI 

Algorithm 3. Let SD of L and the potential Qi{x) for 0 <x< 5 < 7 rbe given. 


1. Construct hi and A by steps 1-3 of Algorithm 1 . 

2. Look for Qi{x) in the form 


Qi{x) 


Qi{x), X <6, 

Qild) + iK,{x — 6)1 x>6 


Then P±{x) = exp {±(Qi(5)(a; — 5) + ^(x — 5)^/)} P±{6), and Hi must satisfy the rela¬ 
tion 

(/ + Hi)-\I - Hi) = P_(7r)i(J - hi){I + hi)-^p-\7r). 

It is equation fl62l) with O in the form 

0{k) = exp{—iK{7r — 6)‘^)l{, U^=U~^. 


Take such k, that det((9(K) -|- /) 7 ^ 0. 

3. Construct Hi via fl63|l . 

4. Put Qo{x) = ho = Hq = 0, X ^ [ 0 , 7 r]. 


Note that if one chooses k bounded (say, |k| < 1), then ||Qi(x)|| < C, where C depends 
only on Qi and does not depend on 6. 

Algorithm 4. Solution of Inverse Problem 1. Let the spectral data of the pencil L be 
given. 

1. Put k := 1, ^0 := 0. 

2. We assume that for x G [0, the potential Qi{x) is already known. Choose the model 
pencil L satisfying fl3T]) . using Algorithm 3. 

3. Step 2 of Algorithm 2. 

4. Choose the largest 5*,, 5k-i < 6k < 7i, such that the main equation flFTD has a unique 
solution for x G {6k-i,6k), hnd this solution. 

5. Steps 4-5 of Algorithm 2 for a: G (5fc_i, 6k). 
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6. Substituting ip{x,p) into ([I]) and ([2]), obtain Qs{x), x G {5k-i-,5k)-, hg, s = 1,2. If 6k = vr, 
construct also Hg, s = 1,2. 

7. If 6k = vr, then terminate the algorithm, otherwise put k := k + 1 and go to step 2. 

Theorem 4. Algorithm 4 terminates after a finite number of steps. 

Proof. Consider the pencil L constructed by Algorithm 3. 

We have P±{x) = P±{x) for x G [0, h]. For x > 6, the following integral representations are 
valid 




Qi{t)P±(t) dt, P±{x) = P±{6)± Qi{x)P±(t) dt. 


Then ||P±(a;)|l, ||P±(a;)|| < C, where the constant C depends only on Qi. Using ([39]), we obtain 


U(a;) = I + Bs{x), \\Bs{x)\\ < Cq^{x - 5). 


Choose an integer s such that vr/s < 1/(2CqJ. Put Xk ■= kir/s, k = 0,s. Then if we 
construct Qi by Algorithm 3 with 6 = Xk-i, we get ||i? 5 (a;)|| < 1/2, x G {xk-i,Xk). Hence 
det r2(x) 7 ^ 0 for X G {xk-i, Xk) and the main equation fl57D is uniquely solvable in this interval. 
It is easy to show that 6k > Xk for all fc = 1, s. Since Xg = vr, the number of steps in 


Algortihm 4 is hnite. 


□ 


6. Multiple poles of the Weyl matrix 

In this section, we consider the general case, when N is not necessary equal zero and the 
Weyl matrix can have multiple poles. For convenience, we assume that for \n\ > N the assertion 
of Theorem [H holds for both collections and 

Let a contour Cq enclose the set {png, pnq}\n\<N and do not contain other eigenvalues pnq, 
pnq, |vr| > N, inside or on the boundary. Then instead of relations flT7|) and (jUj), we get 



m 



|n|>7V q=l 



m 



|n|>A^ q=l 


Substituting ([5|) into int integral in we obtain 



where 



m„ —1 
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n ^ S, k = 0, mn — 1, 

1 Qu—k 

¥^n+fc,l(^) • ^ Pn)i ^n+fc,l(^)P) ^ ^ _ fc)! ^^”-)’ 

n G /c = 0, — 1. 

It is easy to check, that 

Qu Qu+s 

TflCGa ^ r^jCGa ^ 

uniformly with respect to x G [0, tt] and p, 0 on bounded sets, for Sr^x^p) and SRjkix, p,9) 
from the proof of Lemma [SI Hence one can differentiate fIMD by p and 0651) by p and 9. From 
(EH) we get 


d'^ d'^ I d'^ 

- I '^^nfl{x)^Anfl{x, p) - ^(pn,l{x)—An,l{x, p) 
n=l " n=l 


dp^ 


dp'' 


dp'' 


Qv Qv 

EE( 'f{^-)Pnq)TTT'nq^'riqQ^j^D{x,p,Pnq) i Pnq)'^nq ^nq ^D{x, P, Pnq) ^ ■ (66) 


Inl^Af q=l 


‘dp’' 


Relations fHTl) and fH2D can be modihed in a similar way. 
Consider the hnite-dimensional space 

Ho:= {6 = {64f+/:6fceC“}, 


\bo = max _ ||6fc||. 

l<k<K+K 


Denote 

p{x,p)\Bo ■= {pi,o{x),. . . , Pk,o{x), Pl,l{x), . . . , Pj^^^ix)} G Bq. 

Let Gk, k G Z\{0} be the sets |n| > iV, g = l,p, renumerated in a natural way. Instead 
of H, introduce the following Banach space 

= Ho X Hi = {/ = {/fcjfcGZ, /o £ Ho, fk G B{Gk), k G Z\{0}; ||/||_b* := sup ||/fc|| < cxd}. 
Fix X G [0, tt]. Then p(x, p) produces an element of H* by the following rule 


^/>(x) := (p(x,p)|i5, = {'0fc(a;)}fce2, 'ipoi.x) = p{x, p)\Bo, ^/’fc(x) = (p(x, p)b(g,), fc G Z\{0}. 

Similarly ^lJ{x) := <p(x,p)|B,. Then the flSSD and the analogue of flTT]) can be transformed into 
the following equations in the Banach space Bp. 


VL{x)'ip{x) = 'ip{x){I + R{x)), Q{x)iIj{x) = 'ip{x){I — R{x)), 

where I is the identity operator in H*, and R{x), R{x) are linear operators, acting from H* to 

H,. 

According to fISS]) . the operator R{x) acts to an arbitrary element "0 G H in the following 
way: 

'ipR{x) = '^'ipkRkix), 0 = {0fc}fcGZ, (67) 

fcez 


i^oRoix) 


K K 

'IpOjAj^oix, P)-Y '^0,K+jAj,i{x, p) 
j=i i=i 




|B. 


21 

















i)kRk{x) = 


^ ^ '4^k^Pnq)^nq^nqDiS^i P: Pnq) ^ ^ '4^k^Pnq)^riq^nqD{,^: Pi Pnq) 

^iQ'- Pnq^Gk 'R'lQ- Pnq^Gk 




The results of Lemma [HI Theorems [2] and [3] and Remark [Hare valid for this operator as well 
as for one considered in Section 4. Namely, for each hxed x G [0, tt], 


• The series in (I37|) converges absolutely and the operator R{x) is bounded in S*. The proof 
is analogous to Lemma [HI since the inhnite-dimensional part of the considered operators 
is the same (given by Rk,n{,x)). 

• Equation fl37D (considered now as an equation in B^) has a unique solution z{x) G for 
such X G [0, tt] that det ff(a;) 7 ^ 0. To prove this, obtain from fl65l) a relation in similar 
to fl5T|) . Moreover, the operator R{x) is compact in R* and the operator (J + R{x)) is 
invertible on B^. 

• Equation fl57|l in the Banach space R* can be used in Algorithm 4 (with obvious modih- 
cations of the Algorithm) for the recovery of the pencil L from its spectral data SD. 
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